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In the following article we describe our measurement of the proton structure function F 2 and of 
the gluon momentum density xg in ep collisions with the ZEUS detector at HERA in 1993. The 
results for confirm our measurement from the previous year but with much higher statistics 
and show a strong rise towards small x. The gluon momentum density xg is measured for the 
first time in a range from 4 ■ 10~ 4 < x < 10~ 2 at a value of Q 2 = 20 GeV 2 . 



Abstract 



1 Introduction 

In 1993, electrons of 26.7 GeV and protons of 820 GeV were brought into collision at 
HERA. 84 paired bunches out of 210 possible were filled and in addition 10 electron and 6 
proton bunches were left unpaired for background studies. 

The measurements described here were made using the ZEUS detector, which is described 
in detail in [1] and [2]. The main components used for our analyses were the vertex detector 
(VXD), the central tracking detector (CTD), the uranium calorimeter (CAL), a lead scintillator 
beam monitor (C5) and the luminosity detector. 

With the data taken in 1993 a measurement of the proton structure function F% was 
possible with high precision in a region of 7 < Q 2 < 10 4 GeV 2 and x values as low as 3 x 10~ 4 . 
The total integrated luminosity in 1993 was 0.54 pb _1 . 

2 Kinematics 

The processes in deep inelastic scattering (DIS) are described by the following kinematic 
variables. 
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k and k! denote the four-momenta of the incoming and the final state lepton, P is the four- 
momentum of the incoming proton and y is the fractional energy transfer in the proton rest 
frame. 

As the ZEUS detector is a nearly hermetic detector (99.7 %), there exist various methods 
to reconstruct the kinematic variables. For our analysis we have used two different methods: 

1. electron method; (in this method only the scattered electron is observed and its scattering 
angle and energy are measured.) 

2. double angle (DA) method; (this relies on the measurement of two scattering angles, that 
of the scattered electron and that for the final hadronic state. In the naive Quark Parton 
Model the latter is the angle of the scattered quark.) 

Method 1 gives a better resolution in x at low Q 2 while method 2 is less sensitive to 
the absolute calorimeter energy scale and gives better resolution over the whole x. — Q 2 plane. 
Comparing the results obtained by both methods gives a good systematic check of our analysis. 

The range of x and Q 2 that is accessible with these two methods is shown in figure 1. 




Figure 1: Distribution of events in the x and Q 2 plane. Also shown are the bins used and the 
range covered by previous experiments. 
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3 Event Selection 

The event selection is dominated by two selection steps, the trigger and the offline selection. 
The ZEUS trigger is a three step algorithm and has an efficiency greater than 95% for DIS events 
in the x and Q 2 region discussed here. A total of 2 ■ 10 6 events passed the DIS trigger. To obtain 
an event sample with little background and high purity, several cuts were applied. The most 
important are listed below: 

— 35 < 8 < 60 GeV, 6 = 51; Ei(l — cos #.;), where Ei denotes the energy deposit in calorimeter 
cell i and the polar angle of the centre of the cell with respect to the event vertex; (For 
fully contained DIS events 8 m 2E e = 53.4 GeV, whereas photoproduction events peak at 
low values of 8, as the scattered electron remains within the rear beampipe.) 

— Eg > 5 GeV; the energy of the scattered electron should be greater than 5 GeV for a reliable 
electron finding and a low background from photoproduction; 

— llJB > 0.04 to ensure a good measurement of the hadronic angle and of x. 

A total of 46k DIS NC events corresponding to an integrated luminosity of 0.54 pb _1 
survived all cuts. 
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Figure 2: Distribution of the primary measured quantities, a) angle of the scattered electron $' e , 
b) angle of the hadronic final state jjj, c) Q 2 , and d) x. The dots are the data and the shaded 
represent the MC. 

The distributions for the measured and the reconstructed variables for the DA method for 
data and MC are shown in figure 2. The good agreement between the data and MC distributions 
demonstrates our good understanding of the relation between MC and data. The structure 
function used as an input for the MC distribution has been reweighted to the measured structure 
function. 

4 Background 

The main source of background originates from photoproduction events. Figure 3 shows 
the distribution of DIS and photoproduction events in a given x and Q 2 bin after all cuts, except 
for the cut on 8, have been applied. To calculate the number of photoproduction events with 
8 > 35 GeV a fit is made to the data, shown as black dots, in each bin of the x — Q 2 plane. From 
the total number of events in each bin, the number of photoproduction events as calculated 
from the fit, is subtracted. The photoproduction background can rise to 12% at high y, but is a 
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few percent for most of the bins. Other sources of background, like cosmic ray events and halo 
muons, contribute much less and are efficiently removed. 
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Figure 3: Distribution of DIS and photoproduction (php) events in a given x — Q 2 bin. The black 
dots show the data and the solid line represents the fit that is applied to the data. 



5 The Proton Structure Function F 2 

The proton structure function F<i is related to the double differential cross section for 
inclusive ep scattering including photon and Z° exchange via the following formula: 



d 2 a 2ira 



2 



dx dQ 2 xQ 4 



[Y + T 2 (x, Q 2 ) - y 2 F L (x, Q 2 ) + Y-xF 3 (x, Q 2 )] (1 + S r (x, Q 2 )) (2) 



with Y± = 1 =b (1 — y 2 ). Tl is the longitudinal structure function, Tj, is the parity violating 
term, and 6 r is the radiative correction. For comparison with other experiments we measure 
F2 containing only contributions from virtual photon exchange. Therefore we have to absorb 
the contributions from the Z° exchange in a correction factor 8z- This correction is 6% for the 
highest Q 2 bin shown in this analysis and lower for all the other bins. Also the corrections due 
to the longitudinal structure function and the radiative corrections are small (lower than 15% 
and 5% respectively), and can be absorbed in correction factors Sp L , from QCD prescriptions, 
and 6 r , from the HERACLES Monte Carlo. This results in the following equation for F2: 

To derive F2 from the observed double differential cross section, an unfolding method is used. 

The results for the DA method are shown in figure 4. Our most important systematic 
check is the comparison of the values of F2 for the DA and the electron method. These two 
methods use overlapping data samples, but the reconstruction of the kinematic variables is 
largely independent. Within the statistical errors there is very good agreement between the two 
methods. 

6 Gluon Momentum Density 

In zeroth order, F2 is constant in Q 2 and the gluons do not contribute to the cross section. 
The scaling violations, as shown in figure 5, are a higher order process (called leading order, 
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Figure 4: Proton structure function Fi as a function of x. The inner error bars show the statistical 
error only, while the outer bars show the statistical and the systematic error added in quadrature, 
is given in GeV 2 . 



LO, for our purpose). At small values of x, x < 10~ 2 , they are dominated by quark pair produc- 
tion from gluons leading to positive scaling violations with increasing Q 2 . Therefore the gluon 
momentum density, xg(x,Q 2 ), can directly be extracted from the slope dF^j 'dlnQ 2 . 

We compare different methods, the approximate method of Prytz [5], EKL [6], and a fit 
making use of the full NLO GLAP [7] evolution equation. 

The method of Prytz uses the approximation that at low values of x the GLAP evolution 
equation is dominated by the gluonic terms. Consequently the quark splitting functions are 
neglected and only the gluon splitting functions are considered. Using a LO result for P qg , 
the gluon splitting function and making a Taylor expansion of the gluon momentum density 
zg(z, Q 2 ) at z = x/2 results in the following equation in LO: 

(40/27)a< l /47r 

where the gluon density at x is calculated from the slope of F2 at x/2. 

Equation 4 has been extended to NLO, still neglecting quark contributions. In the NLO 
GLAP fit, all the F2 data from ZEUS taken in 1993 have been used. To make a constraint at large 
values of x, NMC data [8] on F2 were included. The fit of the parton momentum distributions 
as a function of x has been performed at a value of Qq = 7 GeV 2 assuming a functional form of 
the parton distribution as used by many parametrisations [4]. 

The power of the rise of the gluon momentum density with decreasing x has been calculated 
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Figure 5: Scaling violation of the proton structure function F2 as a function of Q for fixed 
values of x as given in the figure. The solid lines show a straight line fit, while the dashed lines 
show the result of the global NLO fit. The arrows indicate Q 2 = 20 GeV 2 . 

and found to be xg ~ A ■ x u - JO -o.o4. 

The results from the global fit, from the method of Prytz and from that of EKL are 
shown in figure 6 for Q 2 = 20 GeV 2 . The shaded band shows the error on the global fit. The 
inner error bars on the points show the statistical error only, while the outer error bars show the 
statistical and the systematic error added in quadrature. A good agreement between all methods 
is observed. 



7 Conclusions 

With the 1993 ZEUS data the strong rise of F2 with decreasing x already observed in the 
previous year has been confirmed and is clearly visible over a large range of Q 2 . The measured 
values of F2 are in agreement with the parton distribution parametrisations from GRV, MRSA 
and CTEQ2D'. The observed scaling violations have made it possible to derive the gluon mo- 
mentum density in a new regime of x at Q 2 of 20 GeV 2 . xg(x,Q 2 ) shows a steep rise at low x 
when compared to the measurements from NMC at large x. 
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Figure 6: Gltion momentum density as a function of x at Q 2 = 20 GeV 2 
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